Colloidal glass transition: beyond mode-coupling theory 
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A new theory for dynamics of concentrated colloidal suspensions and the colloidal glass transition 
is proposed. The starting point is the memory function representation of the density correlation 
function. The memory function can be expressed in terms of a time-dependent pair-density corre- 
lation function. An exact, formal equation of motion for this function is derived and a factorization 
approximation is applied to its evolution operator. In this way a closed set of equations for the den- 
sity correlation function and the memory function is obtained. The theory predicts an ergodicity 
breaking transition similar to that predicted by the mode-coupling theory, but at a higher density. 
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There has been a lot of interest in recent years in the 
theoretical description of dynamics of concentrated sus- 
pensions and the colloidal glass transition It has 
been stimulated by ingenious experiments which provide 
detailed information about microscopic dynamics of col- 
loidal particles Q . Due to the abundance of experimen- 
tal data the colloidal glass transition has emerged as a 
favorite, model glass transition to be studied 

One of the conclusions of these studies is the accep- 
tance of the mode-coupling theory (MCT) as the the- 
ory for dynamics of concentrated suspensions and their 
glass transition 0. Historically, this is somewhat sur- 
prising since MCT was first formulated for simple fluids 
with Newtonian dynamics [5j and only afterwards was 
adapted to colloidal systems with stochastic (Brownian) 
dynamics On the other hand, basic approximations 
of MCT are less severe for Brownian systems |jf ■ 

MCT is a theory for correlation functions of slow vari- 
ables, i.e. variables satisfying local conservation laws. 
For Brownian systems there is only one such variable: lo- 
cal density. MCT's starting point is the memory function 
representation of the density correlation function 0, [t| . 
The memory function is expressed in terms of a time- 
dependent pair-density {i.e. four-particle) correlation 
function evolving with so-called projected dynamics. For 
Brownian systems this step is exact [lOj . The central ap- 
proximation of MCT is the factorization approximation 
in which the pair-density correlation function is replaced 
by a product of two time-dependent density correlation 
functions. As a result one obtains a closed, nonlinear 
equation of motion for the density correlation function. 
This equation predicts an ergodicity breaking transition 
that is identified with the colloidal glass transition. MCT 
has also been used to describe, e.g., linear viscoelasticity 
[Tl| . dynamics of sheared suspensions [l^, and colloidal 
gelation [l^. By and large, its predictions agree with 
experimental and simulational results 0, 0] . 

In spite of these successes, MCT's problems are well 
known The most important, fundamental problem is 
that once the factorization approximation is made there 
is no obvious way to extended and/or improve the theory. 



This is most acute for Brownian systems because there 
the density is the only slow mode and thus couplings 
to other modes cannot be invoked! Furthermore, MCT 
systematically overestimates so-called dynamic feedback 
effect. Thus, e.g., it underestimates the glass transition 
volume fraction for a Brownian hard-sphere system (by 
about 10% Q) and overestimates the glass transition 
temperature for a Lennard- Jones mixture (by a factor 
of 2 0]). Finally, MCT cannot describe slow dynamics 
in systems without static correlations [lfij . 

A way to improve upon MCT would be to introduce 
many-particle dynamic variables into the th eory . Such 
an attempt has been made for simple fluids |17| : it was 
argued that these variables (essentially, pair-density fluc- 
tuations) describe clusters of correlated particles. Unfor- 
tunately no quantitative results have been reported based 
on this interesting approach. 

We propose a different way to go beyond MCT. Rather 
than factorizing the pair-density correlation function, we 
derive an exact, formal equation of motion for it [lsj . The 
structure of this equation is very similar to that of the 
equation of motion for the density correlation function; 
"pair" analogues of the usual frequency matrix and the 
irreducible memory function can be identified. The basic 
approximation of our theory is a factorization of the evo- 
lution operator of the pair-density correlation function. 
After this approximation we obtain a closed system of 
equations of motion for the density correlation function 
and the memory function. These equations predict an er- 
godicity breaking transition; for a Brownian hard sphere 
system the glass transition volume fraction, <fi g , is equal 
to .549 (note that 4>f CT = .525, 4>f p « .58). 

Our theory is similar to MCT in that it relies upon an 
uncontrollable factorization approximation. In contrast 
to MCT, it uses this approximation one step later. Thus, 
e.g., our theory preserves the memory function represen- 
tation of the pair-density correlation function while MCT 
approximates the latter by a product of two density cor- 
relation functions. However, as usual in the liquid state 
theory, a priori these features do not guarantee the su- 
periority of our approach as compared to MCT. 
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Our theory starts from the memory function represen- 
tation of the density correlation function, F(k; t), 



where Qi = 1 — P;, and the projection operator Pi reads 



F(k;t) 



1 

N 



(ra(k)exp(fii)n(-k)) 



(1) 



Here N is the number of particles, n(k) is the Fourier 
transform of the density n(k) = J2i e _lk r S and O is the 
iV-particle evolution operator, i.e. the Smoluchowski op- 

erator, = D X), £ • - jSFj) 0, with A) being 
the diffusion coefficient of an isolated Brownian particle, 
(3=1/ (fcfiT), and F; a force acting on particle I. Finally, 
(. . .} denotes the canonical ensemble average; the equi- 
librium distribution stands to the right of the quantity 
being averaged, and all operators act on it as well as on 
everything else. Usually, the memory function represen- 
tation of the Laplace transform of the density correlation 
function, F(k; z), is written as [2(| 



F(k;z) 



S{k) 



(2) 



S{k)(l+M(k-z)) 



where S(k) is the static structure factor and M(k; z) is 
the Laplace transform of the irreducible memory func- 
tion. We re-write @ in a form that will allow us to 
identify the pair analogues of the frequency matrix and 
the memory function. We write a memory function ex- 
pression for the Laplace transform (LT) of F(k; t) 

LT(F(k: t)) = k- (l + M(k;z)0- 1 ) _1 0-k 

x -, , rrr F(k;z). (3) 
(n(k)n(-k)) 

Here 1 denotes a unit 3d tensor, is defined through 
-k • • k = (n(k)Qn(-k)) (note that = 1D N), and 
M(k; z) is the Laplace transform of the current correla- 
tion function evolving with projected dynamics, 



M(k;t) = (j(k)exp(^*)j(-k)) 
where j(k) is a projected current density, 



j(k) =Q n D J2(-^ + l3F 



-ik-Ti 



(4) 



(5) 



q 



><e" 



(8) 



To make connection with the usual form of the mem- 
ory function representation we note that k • • 
k/ (n(k)n(— k)) = Dgk 2 / S(k) is the frequency matrix 
and k • M(k;z)0 -1 • k = M (k; z), where k = k/fc, is the 
irreducible memory function. 

To obtain a convenient expression for M(k;f) in terms 
of a pair-density correlation function we use the following 
exact jjnj equality: 



j(- k ) = E E ^(-ki,-k 2 ) 

(ki,k 2 ) (k 3 ,k 4 ) 

xgQs.i, k 2 ; k 3 , k 4 ) <n 2 (k 3 , k 4 )j(-k)) . (9) 



Here n 2 (ki,k 2 ) is the part of pair-density fluctuations 
orthogonal to the one-particle density fluctuations, 



n 2 (k 1 ,k 2 ) = Q n ^e- ikl ' 



(10) 



Furthermore, in Eq. the sums over k^ < kj+i are 
understood and g denotes the inverse pair-density fluc- 
tuations matrix (it is a pair analogue of 1/ (n(k)n(— k))), 

^ .g(ki,k 2 ;k 3 ,k 4 ) (n 2 (k 3 , k 4 )n 2 (-k 5 , -k 6 )) = 

(k 3 ,k 4 ) 

£ki,k 6 £k 2 ,k 6 - (H) 



Using identity © we can express memory function Q 
in terms of the time-dependent pair-density correlation 
function evolving with one-particle irreducible dynamics, 



In Eq. JSJ Q n = 1 — in, and P n is a projection operator 
on the density subspace, 



4 = ...5>(-q) 



1 



(n(q)n(-q)) 



ra(q) . 



(6) 



Finally, in Eq. ^ zrr 1S the "one-particle irreducible 
Smoluchowski operator" [20| . 

fiirr = ^E^r(^-/3pAQ n , (7) 



F 22 (ki,k 2 ;k 3 ,k 4 ;t) = 

(n 2 (k 1 ,k 2 )exp(O ir '"i)n 2 (-k 3 ,-k4)). (12) 



Rather than factorizing i 22 , we use the projection opera- 
tor method to derive an exact, formal equation of motion 
for this function. The derivation will be given elsewhere 
{jjj; here we present the structure of the final formula 
for the Laplace transform of the time-derivative of the 
pair-density correlation function, P 22 , 
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Lr(F22(k 1 ,k a ;k8,k 4 ;t)) =-111 J2 ] (ki,k 2 )(Z+X(k 5 ,k 6 ;k 7 ,k 8 ;z)O(k 7 ,k 8 ;k 9 ,k 10 )- 1 ) 

i=3 \(k 2i _i,k 2i )/ 



O(k 9 ,kio;kii,ki 2 ) ( ^ J #(ku, k i2 ; k i3 , ki 4 )F 22 (ki 3 , k M ; k 3 , k 4 ; z) (13) 



In Eq. i|13|) ^ denotes a unit 6d tensor, and Ai are 
block matrices, e.g. 



0(k 1 ,k 2 ;k 3 ,k 4 ) = 

Oi 1 (k 1 ,k 2 ;k 3 ,k 4 ) Oi 2 (k 1 ,k 2 ;k 3 ,k 4 ) 
021 (ki , k 2 ; k 3 , k 4 ) 22 (ki , k 2 ; k 3 , k 4 ) 

and the following short-hand notation is used: 



(14) 



(k 1 ,k 2 )0(k 1) k 2 ;q 1 ,q 2 ) 



qi 
q2 



Y^k,-0, r q r (15) 



1,3 



O and Ai are the pair analogues of and M (compare 
Eqs. (|2J and (fT3l '); in particular 



(k^k^kx^k^k. 



^(ki.kajn*" n 2 (-k 3 ,-k 4 )) 



(16) 



and M.ij are pair-current correlations evolving with a 
two-particle irreducible evolution operator £l 2lrr , e.g., 

Xn(ki,k 2 ;k 3 ,k 4 ;t) = 
<j 2 (k 1 ,k 2 )exp(0 2 ^)j 2 (-k 3 ,-k 4 )) (17) 



where, e.g. 



j2(-k 3 ,-k 4 ))=g„ 2 I)o^^Q 1 e ik 3 .r !+ ik 4 .r m ) 



(18) 

Explicit formulae for O and Ai (including definitions of 
fl 2trr and Q n2 ) will be given elsewhere . 

The main approximation of our theory is factorization 
of the evolution operator for i 7 ' 22 . Within this approxi- 
mation the diagonal blocks of O and Ai are given by 

0n(ki,k 2 ;k 3 ,k 4 ) = 22 (k 2 ,ki;k4,k 3 ) - 
NOS(k 2 )S klM d k2 ,i H , (19) 
Xn(ki,k 2 ;k 3 ,k 4 ;t) = X 22 (k 2 , ki; k 4 , k 3 ; t) = 
^M(k i; t)F(fc 2 ;t)5 klik3 5 k:! , k4 , (20) 

and the off-diagonal ones vanish. Consistently, we also 
factorize g and F 22 (t = 0). 

Using (|19H2U[) we can express F 22 in terms of the den- 
sity correlation function and the memory function (note 
that F 22 does not factorize for t > 0). Substituting F 22 
into the formula for the memory function and using con- 
volution approximation for static vertices 0, 0| we get 



M(k;z) 



nDo f dkidk 2 



<5(k-ki-k 2 ) k- (c(fci)ki +c(fc 2 )k 2 



S(k 1 )S(k 2 ) 



z + 



l+LT(M(fe 1 :t)_F(fc 2 ;i)/S(fc 2 )) 



+ (!<-> 2) 



, (21) 



where n is the density and c(k) is the direct correlation 
function. Eqs. J5J) and l|21[) determine time dependence 
of density correlations and the memory function. 

Eqs. @ and l|21[l predict an ergodicity breaking tran- 
sition. In the non-ergodic regime F(k;t) has a non- 
zero long-time limit, lim^oo F(k; t) — f(k)S(k), where 
f(k) is called a non-ergodicity parameter. It follows 



from Eq. J5J that in this regime also the memory func- 
tion has a non-zero long-time limit, lim^oo M(k; t) — 
m(k)Dok 2 1 S(k), and that f(k) and m{k) are related by 



/(*) 



m{k). 



(22) 



1 - /(*) 

Using I|21I22II we get a self-consistent equation for f(k): 



f(k) = jn 

l-f(k) 2P J (27T) 



rfkidk 2 

3-d(k - ki - k 2 ) 



k- (c(fci)ki +c(/c 2 )k 2 



r 



S{k)S(k 1 )S(k 2 )f{k 1 )f{k 2 ) 
l + (l-/(fcx))(l-/(fc 2 ))- 



(23) 



One should note that the right-hand-side of an analogous self-consistent equation derived from MCT has a similar 
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FIG. 1: Non-ergodicity parameter. Lines: theoretical pre- 
dictions at the ergodicity breaking transition; solid line: our 
theory, <j> g = .549; dashed line: MCT, 4>f CT = .525. Sym- 
bols: experimental data taken at <j> = .563 [l4ll . 

form; the difference is that within MCT the right-hand- 
side is a quadratic functional of f(k) |5j whereas in the 
present approach it includes terms of all orders in f(k). 

For low enough densities Eq. (|2*5)l has only trivial so- 
lutions (i.e. f(k) = 0). For the hard-sphere interaction 
a non-trivial solution appears at n 9 7rer 3 /6 = <f> g = .549. 
Qualitatively, the ergodicity breaking transition is simi- 
lar to that predicted by MCT: f(k) has a jump at the 
transition. Also, f(k) at the transition is similar to that 
of MCT at the MCT transition, $\ lCT = .525 (Fig. 1). 

The factorization approximation proposed here is the 
simplest possible one. There are two ways to improve 
upon it. First, one could try to include in an approxi- 
mate way the off-diagonal blocks of A4. To this end one 
could express them in terms of a triple-density correlation 
function and then factorize this function into a product of 
three density correlation functions. Second, since the fre- 
quency matrix involves only static correlations, one could 
try to include it in a more sophisticated way. For exam- 
ple, one could include two-particle dynamics exactly |l8| . 
The second extension could describe glassy dynamics in 
systems without static correlations [l6|. 

To summarize, we proposed a new theory for dynamics 
of concentrated suspensions and the colloidal glass tran- 
sition. The theory goes beyond MCT in that includes, in 
an approximate way, time-dependent pair-density fluctu- 
ations. In contrast to an earlier approach [lj], the present 
one uses pair-density correlation function evolving with 
one-particle irreducible dynamics. The new theory pre- 
dicts an ergodicity breaking transition similar to that of 
MCT, but at a higher density. 

The author benefited from inspiring discussions with 
Rolf Schilling and Arun Yethiraj; support by NSF Grant 
No. CHE-0111152 is gratefully acknowledged. 
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